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Abstract. Starting from an effective Skyrme interaction we study effects of phonon-phonon coupling on the 
low-energy electric dipole response in 124 ~ 132 Sn. The QRPA calculations are performed within a finite rank sepa- 
rable approximation. The inclusion of two-phonon configurations gives a considerable contribution to low-lying 
strength. Comparison with available experimental data shows a reasonable agreement for the low-energy El 
strength distribution. 



1 Introduction 



Exotic nuclear collective excitations, like the pygmy dipole 
resonance, represent a subject of intense investigations du- 
ring last decades, see for example B1I2I3I41 . The structure 
and dynamics of low-energy dipole strength, also referred 
to as pygmy dipole resonance (PDR), has extensively been 
investigated using a variety of theoretical approaches and 
models |5|. Recent studies have made use of the Hartree- 
Fock (HF) plus random phase approximation (RPA) [6], 
the Hartree-Fock-Bogoliubov (HFB) model plus quasipar- 
ticle RPA (QRPA) [7 8|, QRPA plus phonon coupling g), 
the quasiparticle phonon model (QPM) including complex 
configurations I1I10I . the relativistic RPA [11] and QRPA 
[12|. Also the quasiparticle time blocking approximation 
(QTBA) has been used either in a non-relativistic frame- 
work [8] or with relativistic Lagrangians (RQTBA) H 1 31 1411 . 

One of the successful tools for describing the PDR is 
the QRPA with the self-consistent mean-field derived from 
Skyrme effective nucleon-nucleon interactions. Such an ap- 
proach describes the properties of the low-lying states less 
accurately than more phenomenological ones, but the re- 
sults are in a reasonable agreement with experimental data 
(Ref. [5] and references therein). Due to the anharmonic- 
ity of vibrations there is a coupling between one-phonon 
and more complex states |15 16]. The main difficulty is 
that the complexity of calculations beyond standard QRPA 
increases rapidly with the size of the configuration space, 
so one has to work within limited spaces. Using a finite 
rank separable approximation suggested in 11711811911 for 
the residual interaction resulting from Skyrme forces one 
can overcome this difficulty. In this paper we study the 
properties of the low-lying and high-lying electric dipole 
strength in the neutron-rich Sn isotopes. The couplings be- 
tween one- and two-phonon components in the wave func- 
tions of excited states are taken into account. 
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2 The method 

The calculations are performed by using the SLy4 11201 in- 
teraction in the particle-hole (p-h) channel and a density- 
dependent zero-range interaction in the particle-particle (p- 
p) channel. Spherical symmetry is assumed for the HF 
ground states. The strength of the surface-peaked zero-ran- 
ge pairing force is taken equal to -940 MeV fm 3 in con- 
nection with the soft cutoff at 10 MeV above the Fermi 
energy as introduced in Ref. |18|. This value of the pair- 
ing strength is fitted to reproduce the experimental pairing 
energies for both protons and neutrons. 

The residual interaction in the p-h channel and in 
the p-p channel can be obtained as the second deriva- 
tive of the energy density functional with respect to the 
particle density p and the pair density p, respectively. Fol- 
lowing Ref. ifTTl we simplify by approximating it by 
its Landau-Migdal form. Moreover we neglect the I - 1 
Landau parameters (Landau parameters with / > 1 are 
equal to zero in the case of Skyrme interactions). In this 
work we study only normal parity states and one can ne- 
glect the spin-spin terms since they play a minor role. The 
two-body Coulomb and spin-orbit residual interactions are 
also dropped. The expressions for F^, F ph and Fq P , Fq P 
can be found in Ref. ED and in Ref. IU8l . respectively. 
The Landau parameters Fq, Go, FLG q expressed in terms 
of the Skyrme force parameters II2TI depend on kp. As is 
pointed out in our previous works 11 1 71 1 81 one needs to 
adopt some effective value for kp to give an accurate rep- 
resentation of the original p-h Skyrme interaction. For the 
present calculations we use the nuclear matter value for kp. 

In describing the giant dipole resonance (GDR) we sho- 
uld exclude the spurious state due to the center-of-mass 
motion of the nucleus. The spurious state admixture can 
be present as a component in each of the wave functions of 
the excited states. The spurious 1~ state is excluded from 
the excitation spectra by introducing the effective neutron 
q n = -Z/A and proton q p = N/A charges 1151 . An alterna- 
tive way of eliminating the spurious state is to orthogonali- 
ze it with respect to all physical states [22 1 . At was shown 
|6| that eliminating the spurious state by means of effective 
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Table 1. Values of the centroid energy and width calculated 
within the QRPA or QRPA plus phonon-phonon coupling (2PH) 
in comparison with the experimental values (Expt.) taken from 
Ref . 1 2 1 . The chosen energy interval is 1 1 -20 Me V. 



Nucleus 


E (MeV) 






r (MeV) 






QRPA 


2PH 


Expt. 


QRPA 


2PH 


Expt. 


I24 Sn 


16.4 


16.3 


15.3 


4.4 


4.7 


4.8 


126 Sn 


16.2 


16.2 




4.4 


4.7 




128 Sn 


16.1 


16.0 




4.7 


4.7 




130 Sn 


15.8 


15.7 


15.9(5) 


4.8 


4.8 


4.8(1.7) 


132 Sn 


15.5 


15.4 


16.1(7) 


4.9 


5.0 


4.7(2.1) 



charges or orthogonalizing it to all physical states lead to 
very similar results. 

We take into account the coupling between the one- 
and two-phonon components in the wave functions of ex- 
cited states. Thus, in the simplest case one can write the 
wave functions of excited states as 1161 



W V {JM) 



V i 



AihA-%12 



10). (1) 



where | 0) is the phonon vacuum, Q*. t | 0) is the phonon 
creation operator and v labels the excited states. The coeffi- 
cients Ri(Jv), P^liJv) and energies of the excited states E v 
are determined by solving the corresponding secular equa- 
tion of Ref . 0U . 

The two-phonon configurations of the wave function 
([]]) are constructed from natural-parity phonons with mul- 
tipolarities A = 1 , 2, 3, 4, 5. All dipole excitations with ener- 
gies below 35 MeV and 15 most collective phonons of 
A = 2,3,4,5 multipolarity are included in the wave func- 
tion dTJ. It is found that extending the model space for one- 
phonon configurations does not change much the calcu- 
lated energies and transition probabilities. 



3 Results and discussion 

In Figure Q] the calculated dipole spectra for 130 Sn are 
shown. The right part of the figure shows the dipole strength 
function up to 26 MeV. The left panel shows the low-lying 
parts of the corresponding spectrum below 12 MeV. Figu- 
re Q] (a) displays the experimental B(El) distribution. The 
black circles with error bars are the experimental data. In 
order to quantify this resonance-like structure the data are 
fitted with a Lorentzian distribution [ 2 3 1 . Calculations wit- 
hin QRPA are shown in Fig. Q] (b), and the QRPA plus 
phonon-phonon coupling (2PH) results are presented in 
Fig. Q] (c). In the figure, the calculated El strength dis- 
tributions are folded out with a Lorentzian distribution of 
1 MeV width. The general shapes of the GDR obtained in 
the 2PH are rather close to those observed in experiments. 
This demonstrates the quality of a description within 2PH 
in comparison with QRPA. We conclude that the main mech- 
anisms of the GDR formation in I30 Sn are taken into ac- 
count correctly and consistently in the 2PH approach. 

Let us discuss results for the low-energy El strength in 
more detail (see in the left part of Fig . |T} . For the low-lying 




Fig. 1. B(El) strength distributions in 130 Sn. (a) Experimental 
data 1 2 3 1; (b) QRPA results; (c) QRPA plus phonon-phonon cou- 
pling results. 



Table 2. Mean energies E and summed B(El) values for the low- 
energy dipole states in the excitation energy range below 1 1 MeV. 
The experimental values (Expt.) are taken from Refs. [2 3 4|. To 
compare with experimental data in 124 Sn we choose the energy 
interval £ < 10 MeV (4). 



Nucleus 


E (MeV) 






£ B(El) (e 2 fm 2 ) 




QRPA 


2PH 


Expt. 


QRPA 


2PH 


Expt. 


124 Sn 


9.7 


9.1 


6.97 


0.86 


0.59 


0.379(45) 


126 Sn 


10.1 


10.0 




1.82 


1.86 




128 Sn 


10.0 


10.0 




1.63 


1.78 




130 Sn 


10.0 


10.0 


10.1(7) 


1.40 


1.80 


2.4(7) 


132 Sn 


9.9 


9.9 


9.8(7) 


1.27 


1.42 


1.3(8) 



part of the dipole spectrum, the QRPA calculation predicts 
two rather pronounced states around 9.5 MeV. The neu- 
tron transition densities of these levels are dominated out- 
side the nuclear surface. That corresponds to the vibrations 
of a neutron skin against a proton-neutron core. One can 
see from Fig. Q] (b) that the next fairly collective state is 
around 1 1.6 MeV. Using the analysis of the transition den- 
sities we found that states above this energy have an inter- 
mediate behavior. Increasing further the excitation energy 
we observe the low-energy GDR tail. The corresponding 
2PH strength distributions in Fig.[T](c) show many states 
with comparable strength in the energy region below 1 1 
MeV. The transition densities of these states, originating 
the fragmentation of the QRPA pygmy mode, have a be- 
havior which is very similar to the behavior of the initial 
QRPA state. However, we can state that the low-lying tran- 
sition densities are also dominated by the neutron contribu- 
tion outside the nuclear surface. Strictly speaking, a correct 
comparison of the integrated pygmy dipole strength with 
experimental data in this region is still an open problem, 
because of loss of experimental sensitivity 112131 . 

Now we consider the integral characteristics of the 
GDR. The GDR calculated characteristics are displayed in 
TableQ] They are compared with the experimental data (2). 
The centroid energy E and the spreading width F are given 
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by the following expressions: 



mo 



i mo 



(2) 



where m,t = 2 B(El)E k are the energy-weighted moments. 
The integral characteristics of £1 strength function in these 
nuclei have been calculated for the energy interval 11-20 
MeV, which exhausts between 100 and 110% of the Tho- 
mas-Reiche-Kuhn energy-weighted sum rules (14.8NZ/A 
MeV e 2 fm 2 ). 

First of all we checked our approach for 124 Sn and ob- 
tained (see Table Q} a reasonable agreement with avail- 
able experimental data 0. In 124 Sn we find that E = 16.4 
MeV for the QRPA and E = 16.3 MeV for the 2PH, even 
though the average energy of the giant dipole resonance 
is by 1 MeV higher than its experimental value. This is 
in agreement with the RRPA calculations [11] where the 
centroid energy for GDR is located at 15.86 MeV in the 
interval 10-25 MeV. For comparison, the QTBA 00 and 
RQTBA Ql calculations give 15.3 MeV and 15.05 MeV, 
respectively. We note that in the QTBA the 0-30 MeV 
summation interval is used while in the RQTBA this inter- 
val is 10 - 22.5 MeV. One can see from Table [TJ that the 
inclusion of the two-phonon terms results in an increase of 
the resonance width from 4.4 to 4.7 MeV. This value is in 
agreement with the QTBA calculations (4.6 MeV) 0. 

The QRPA results for the neutron-rich 130 Sn are pre- 
sented in Fig.Q](b). The centroid value of the GDR is 15.8 
MeV. The results of our calculations for the El resonance 
taking into account the two-phonon terms are shown in 
Fig. Q](c). Here, the GDR centroid energy is 15.7 MeV, 
that is in agreement with another theoretical work using 
the RRPA ifTfl (E = 15.78 MeV in the 10-25 MeV in- 
terval). The RQRPA and RQTBA results are very similar 
(E = 15.13 MeV for RQRPA and E = 14.66 MeV for 
RQTBA) in the interval 10-22.5 MeV fl2). For the GDR 
width we obtain r = 4.8 MeV that can be compared with 
the other theoretical values 3.49 MeV in the RQRPA and 
4.74 MeV in the RQTBA. The experimental value of the 
width is equal to 4.8(1.7) MeV 00- 

For 132 Sn, values calculated by QRPA and 2PH for the 
GDR energy centroid are E = 15.5 MeV and E = 15.4 
MeV, respectively. The experimental value of the GDR 
width for 132 Sn is 4.7(2. 1) MeV [2| and this is in agreement 
with our results. The inclusion of the two-phonon terms 
gives a small increase of the resonance width from 4.9 to 
5.0 MeV. Our results for 132 Sn are qualitatively similar to 
calculations in Ref. [9] where the Skyrme interaction SIII 
has been used. The peak energy and width with the phonon 
coupling effects are 15.5 MeV and 5.8 MeV, respectively. 
It is worth to note the GDR centroid energy in the inter- 
val 8-25 MeV is 14.7 MeV in the RRPA and it is 14.4 
MeV in the RRPA with the inclusion of particle-phonon 
coupling [12]. The inclusion of particle-phonon coupling 
in the RRPA calculation results in an increase of the reso- 
nance width from 3.3 to 4.0 MeV. 

One can see that the agreement with the experimen- 
tal data for 124 Sn and 130,132 Sn is good. This gives us con- 
fidence to use the model for a prediction of the integral 
characteristics of the GDR in 126 ' 128 Sn. 

In the case 126 Sn, the centroid energy for the 2PH case 
has the same value as for the QRPA calculation 16.2 MeV. 



One can see that for 128 Sn the centroid energy with 2PH 
is 16.0 MeV, to be compared with the value 16.1 MeV of 
QRPA. The values of energy centroids for 126 Sn and 128 Sn 
are rather close to the empirical systematics 31.2A -1 ^ 3 + 
20.6A- 1/6 [231. The empirical predictions 120A" 2/3 (24) 
for the GDR widths are reproduced well in both cases. 

As it is seen from the Table [TJ our calculations show a 
slight lowering of the GDR peak energy and an increase of 
the GDR width when passing from 124 Sn to 132 Sn. These 
results are in an agreement with the calculations of Litvi- 
nova 11311411 and Piekarewicz [11]. The role of anharmonic 
effects, which mainly determine the value of the GDR 
width, is decreasing when passing from non-magic to magic 
nuclei. 

Finally, let us discuss the low-energy El strength. In 
contrast to the GDR, the effects of phonon-phonon cou- 
pling is supposed to be more serious. The integral charac- 
teristics of the strength distribution of low-energy dipole 
states before 1 1 MeV are given in Table [2] Note that to 
compare with experimental data in 124 Sn we choose the 
energy interval E < 10 MeV (4). 

In the case of 124 Sn, the QRPA mean energy is equal 
to 9.7 MeV and taking into account phonon-phonon cou- 
pling gives rise to a decrease of the PDR energy by 0.6 
MeV, while experiment gives E — 6.97 MeV |@). For com- 
parison, the QTBA QD and RQTBA ffTU calculations give 
8.7 MeV and 8.15 MeV, respectively. The experimental 
data for the 4-10 MeV interval give for the integrated PDR 
strength a value %B(E1) = 0.379(45) e 2 fm 2 , while our 
calculations give values 0.86 and 0.59 e 2 fm 2 within the 
QRPA and 2PH, respectively. The calculated total QPM 
ifTOl dipole strength in the PDR energy range E = 5.7-7.2 
MeV is 0.324 e 2 fm 2 . This is rather close to the experimen- 
tally deduced strength. On the other hand, the QTBA 1 8 1 
and RQTBA lfT4l results are very similar for the interval 
below 10 MeV (£fi(£l) = 3.0 e 2 fm 2 within the QTBA 
and 2 B(El) = 3.2 e 2 fm 2 within the RQTBA). 

Values predicted by QRPA and 2PH for the PDR en- 
ergy centroid in 126 Sn are E = 10.1 MeV and E = 10.0 
MeV, respectively. The inclusion of the two-phonon terms 
results in a slight increase of the summed B(El) from 1.82 
to 1.86 e 2 fm 2 . In the case of 128 Sn, the calculated summed 
B(El) strength in the energy range below 1 1 MeV amounts 
to 1.63 e 2 fm 2 (in the QRPA) and 1.78 e 2 fm 2 (in the 2PH). 
The centroid energy is 10.0 MeV in both cases. 

Our calculations give for 130 Sn a total dipole strength 
2 B(El) = 1.40 e 2 fm 2 for the QRPA and £ B(El) = 1.80 
e 2 fm 2 for the 2PH. The summation is performed for the 
dipole states below 11 MeV. The experimental value is 
2 B(El) = 2.4(7) e 2 fm 2 0. On the other hand, the cor- 
responding mean energy of the PDR is 10 MeV (it is re- 
produced in both cases). It may be compared with the ex- 
perimental value 10.1(7) MeV. The RQRPA and RQTBA 
results in the interval below 10 MeV lTT3l are very simi- 
lar £ B(El) = 4.04 e 2 fm 2 (RQRPA) and £ fl(£l) = 3.44 
e 2 fm 2 (RQTBA). The PDR centroid in the RQTBA is equal 
to 7.5 MeV. For comparison, the RRPA calculations give 
E = 7.91 MeV in the same interval ifTTIl . As one can see 
from Figure QJ the inclusion of the two-phonon terms re- 
sults in an essential increase of the PDR width - rpoR- For 
the QRPA case we get rpaR =1.0 MeV that can be com- 
pared with the value rp^R = 1.8 MeV when one takes 
into account the phonon-phonon coupling. Only an upper 
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experimental limit for the PDR width could be deduced 
because of the finite energy resolution in experiment on 
Coulomb dissociation (2). This upper limit for the PDR 
widths is r PDR < 3.4 MeV. 

Now we discuss the low-lying dipole response in 132 Sn 
for the energy interval below 1 1 MeV. In RPA calculations 
the mean energy is 9.9 MeV. Taking into account phonon- 
phonon coupling gives the same value for the PDR energy, 
while experiment gives E = 9.8(7) MeV [2|. We notice 
that the inclusion of phonon-phonon coupling in our cal- 
culations leads to an increase of the 2 B(E1) values from 
1.27 to 1.42 e 2 fm 2 . In the QTBA calculations [8) the cen- 
troid energy of PDR is located at 8.9 MeV for the energy 
interval below 10 MeV. This value is by 0.9 MeV lower 
than its experimental one. For comparison, the RQTBA 
1 14| calculations give 7.3 MeV for the same energy inter- 
val. Our calculation shows that the inclusion of the two- 
phonon terms results in an increase of the resonance width 
from 1 .2 to 2.0 MeV. For the widths of the PDR an upper 
limit is r PDR < 2.5 MeV J2). 

Table [2] demonstrates that there is a distinct difference 
between characteristics of the low-lying strength integral 
in 124 - 132 Sn. Namely, the phonon-phonon coupling con- 
tribution to 2 B{El) is small for nuclei in 126,128,130,132^ 
while it is rather important in the 124 Sn. Moreover, for nu- 
cleus such as 124 Sn the PDR is almost completely defined 
by complex configurations in the energy interval below 10 
MeV. 

The QRPA predicts a monotonic decrease of the PDR 
centroid with mass number. The value of the total PDR 
strength decreases in going to the heavier tin isotopes, up 
to 132 Sn. Other QRPA calculations with Skyrme forces Q 
give similar tendencies for the summed B(El) values but 
our YjB(EI) values are somewhat smaller. The influence 
of coupling between one- and two- phonon terms in the 
wave functions ([1) only leads to the fragmentation of El 
strength. On the other hand, as one can see from Figure 
Q] the inclusion of the two-phonon terms results in an es- 
sential increasing of the PDR width. Also from the Table 
|2] one can see that dipole strength is transferred from the 
GDR region to the low-energy region. In summary, our E 
and 2 B(El) values are in good agreement with the exper- 
imental values, except for 124 Sn. 



4 Summary 

Starting from the Skyrme mean-field calculations, the prop- 
erties of the electric dipole strength in tin isotopes are stud- 
ied by taking into account the coupling between one- and 
two-phonons terms in the wave functions of excited states. 
The finite rank separable approach for the QRPA calcula- 
tions enables one to reduce considerably the dimensions 
of the matrices that must be inverted to perform nuclear 
structure calculations in very large configuration spaces. 

Neutron excess effects on the PDR excitation energy 
and transition strength have been investigated for I24 ~ I32 Sn. 
We find that the centroid energy and summed strength de- 
creases when increasing the neutron number. The inclusion 
of the two-phonon configurations results in an increase of 
the widths of the PDR and GDR. This is in agreement with 
available experimental data. For i 26 -i 28 Sn the properties of 
the PDR are predicted. 
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